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Abstract. We compare the option pricing formulas of Louis Bachelier and 
Black-Merton-Sclioles and observe - tiieoretically as well as for Bachelier's 
original data - that the prices coincide very well. We illustrate Louis Bache- 
lier's efforts to obtain applicable formulas for option pricing in pre-computer 
time. Furthermore we explain - by simple methods from chaos expansion - 
why Bachelier's model yields good short-time approximations of prices and 
volatilities. 



1. Introduction 

It is the pride of Mathematical Finance that L. Bacheher was the first to analyze 
Brownian motion mathematically, and that he did so in order to develop a theory 
of option pricing (see 2J). In the present note we shall review some of the results 
from his thesis as well as from his later textbook on probability theory (see [31), 
and we shall work on the remarkable closeness of prices in the Bachelier and Black- 
Merton-Scholes model. 

The "fundamental principle" underlying Bachelier's approach to option pricing 
is crystallized in his famous dictum (see [2], P-34) 

"L'esperance mathematique du speculateur est nul", 

i.e. "the mathematical expectation of a speculator is zero" . His argument in favor 
of this principle is based on equilibrium considerations (see [2] and similar to 

what in today's terminology is called the "efficient market hypothesis" (see [10]), 
i.e. the use of martingales to describe stochastic time evolutions of price movements 
in ideal markets. L. Bachelier writes on this topic (see the original french version 
in [2], p. 31). 

"It seems that the market, the aggregate of speculators, can believe 
in neither a market rise nor a market fall, since, for each 
quoted price, there are as many buyers as sellers." 

The reader familiar with today's approach to option pricing might wonder where 
the concepts of "risk free interest rate" and "risk neutral measure" have disap- 
peared to, which seem crucial in the modern approach of pricing by no arbitrage 
arguments (recall that the discounted price process should be a martingale under 
the risk neutral measure). As regards the first issue L. Bachelier applied his "fun- 
damental principle" in terms of "true" prices (this is terminology from 1900 which 
corresponds to the concept of forward prices in modern terminology), since all the 
payments involved (including the premium of the option) were done only at matu- 
rity of the contracts. See [2] for an explicit description of the trading rules at the 
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bourse de Paris in 1900. It is well-known that the passage to forward prices makes 
the riskless interest rate disappear: in the context of the Black-Merton-Scholes for- 
mula, this is what amounts to the so-called Black's formula (see ^y)- As regards the 
second issue L. Bachelier apparently believed in the martingale measure as the his- 
torical measure, i.e. for him the risk neutral measure conincides with the historical 
measure. For a discussion of this issue compare [lOj . 

Summing up: Bachelier's "fundamental principle" yields exactly the same recipe 
for option pricing as we use today (for more details we refer to the first section 
of the St. Flour summer school lecture |llj): using forward prices ("true prices" 
in the terminology of 1900) one obtains the prices of options (or of more general 
derivatives of European style) by taking expectations. The expectation pertains 
to a probability measure under which the price process of the underlying security 
(given as forward prices) satisfies the fundamental principle, i.e. is a martingale in 
modern terminology. 

It is important to emphasize that, although the recipes for obtaining option 
prices are the same for Bachelier's as for the modern approach, the arguments 
in favour of them are very different: an equilibrium argument in Bachelier's case 
as opposed to the no arbitrage arguments in the Black-Merton-Scholes approach. 
With all admiration for Bachelier's work, the development of a theory of hedging 
and replication by dynamic strategies, which is the crucial ingredient of the Black- 
Merton-Scholes-approach, was far out of his reach (compare [IT] and section 2.1 
below) . 

In order to obtain option prices one has to specify the underlying model. We fix 
a time horizon T > 0. As is well-known, Bachelier proposed to use (properly scaled) 
Brownian motion as a model for forward stock prices. In modern terminology this 
amounts to 

(1.1) ■.^So + a^'Wt, 

for < t < r, where {Wt)o<t<T denotes standard Brownian motion and the super- 
script B stands for Bachelier. The parameter > denotes the volatility in the 
Bachelier model. Notice that in contrast to today's standard Bachelier measured 
volatility in absolute terms. In fact, Bachelier used the normalization H = 
and called this quantity the "coefficient of instability" or of "nervousness" of the 
security S. The reason for the normalisation H — is that H^/T then equals 
the price of an at the money option in Bachelier's model (see [2]) 

The Black-Merton-Scholes model (under the risk-neutral measure) for the price 
process is, of course, given by 

(1.2) sr-SoeM^'^'W^^^^^t), 

for < t < T. Here a^^ denotes the usual volatility in the Black-Merton-Scholes 
model. 

This model was proposed by P. Samuelson in 1965, after he had - led by an 
inquiry of J. Savage for the treatise [3] - personally rediscovered the virtually for- 
gotten Bachelier thesis in the library of Harvard University. The difference between 
the two models is somewhat analogous to the difference between linear and com- 
pound interest, as becomes apparent when looking at the associated Ito stochastic 
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differential equation, 

This analogy makes us expect that, in the short run, both models should yield 
similar results while, in the long run, the difference should be spectacular. Fortu- 
nately, options usually have a relatively short time to maturity (the options con- 
sidered by Bachelier had a time to maturity of less than 2 months) . 

2. Bachelier versus Black-Merton-Scholes 

We now have assembled all the ingredients to recall the derivation of the price of 
an option in Bachelier's framework. Fix a strike price K, a horizon T and consider 
the European call C, whose pay-off at time T is modeled by the random variable 

Applying Bachelier's "fundamental principle" and using that S^^ is normally 
distributed with mean 5*0 and variance {a^Y'T, we obtain for the price of the 
option at time t = Q 

^ E[{S^ - K)+] 

f°° 1 x'^ 

<") - J.J'' + - - ''■';;5^ 

(2.2, .(S.^Kw'^)^a->Vf,(,^), 

where 0(a;) = exp(— ^) denotes the density of the standard normal distribu- 
tion. We applied the relation (j>'{x) = —x(f>{x) to pass from (|2.ip to (|2.2p . For 
details see, e.g. , [5|. 

For further use we shall need the very well-known Black-Merton-Scholes price, 
too, 

tBS\2 



(2.3) = / {Socxpi-^-^-^ + a'^'^Vrx) - K)+-=exp{--)dx 

J-oo ^ v27r ^ 

/•oo ( BS\2rp -, 2 

(2.4) = y.„, K.,(.^^,^. ^SoeM~-^^ + 'y^'VTx)^K)^cxp{--)dx 



loe ^ + li(jBS\2rp ^ So _ It BS\2rp 

(2-5) -^o$( ^ )-K<^ r^ J" ^ )■ 



Interestingly, Bachelier explicitly wrote down formula (j2.1[) . but did not bother 
to spell out formula ()2.2p . see [H p. 50]. The main reason seems to be that at 
his time option prices - at least in Paris - were quoted the other way around: 
while today the strike prices K is fixed and the option price fluctuates according 
to supply and demand, at Bachelier's times the option prices were fixed (at 10, 20 
and 50 Centimes for a "rente", i.e., a perpetual bond with par value of 100 Francs) 
and therefore the strike prices K fluctuated. What Bachelier really needed was the 
inverse version of the above relation between the option price Cq and the strike 
price K. 
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Apparently there is no simple "formula" to express this inverse relationship. 
This is somewhat analogous to the situation in the Black-Merton-Scholes model, 
where there is also no "formula" for the inverse problem of calculating the implied 
volatility as a function of the given option price. 

We shall see below that L. Bachelier designed a clever series expansion for C(f 
as a function of the strike price K in order to derive (very) easy formulae which 
approximate this inverse relation and which were well suited to pre-computer tech- 
nology. 

2.1. At the money options. Bachelier first specializes to the case of at the money 
options (called "simple options" in the terminology of 1900), when 5*0 = K. In this 
case (|2.2p reduces to the simple and beautiful relation 

As explicitly noticed by Bachelier, this formula can also be used, for a given price 
C — Cq of an at the money option with maturity T, to determine the "coefficient 
of nervousness of the security" H = -^j^, i-e., to determine the implied volatility 
in modern language. Indeed, it suffices to normalize the price Cq by ^/T to obtain 
H = We summarize this fact in the subsequent proposition. For convenience 
we phrase it rather in terms of than of H . 

Proposition 1. The volatility in the Bachelier model is determined by the price 
Cq of an at the money option with maturity T by the relation 

(2.6) - 




In the subsequent Proposition, we compare the price of an at the money call 
option as obtained from the Black-Merton-Scholes and Bachelier's formula respec- 
tively. In order to relate optimally (see the last section) Cq and Cq^ we choose 
= Sq(j and a^^ — a for some constant ct > 0. We also compare the implied 
volatilities, for given price Cq of an at the money call with maturity T, in the Bache- 
lier and Black-Merton-Scholes model. We denote the respective implied volatilities 
by cr^ and a^^ and discover that the implied Bachelier volatility estimates the 
Black-Scholes implied volatility quite well at the money. 

Proposition 2. Fix ct > 0, T > and Sq — K (at the money), let and 
= 5*0(7 and denote by C^ and C^^ the corresponding prices for a European call 
option in the Bachelier il.l]) and Black-Merton-Scholes model il.'^) respectively. 
Then 

(2.7) < - Ci' < -^a'Ti = OiiaVrn 

The relative error can be estimated by 

r<B r<BS rp 

(2.8) ^° < 

Conversely, fix the price < Cq < Sq of an at the money option and denote by 
the implied Bachelier volatility and by a^^ the implied Black-Merton-Scholes 
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volatility, then 
(2.9) 



Proof, (compare [S| and [TT]). For Sq = K , we obtain in the Bachelier and Black- 
Merton-Scholes model the foUowing prices, respectively, 



/2n 



Hence 



BS 




{1 - exp{^ — ))dy 



Vt 



So 



OiiaVTy) 



since e 

^0 



> I + y for all y, so that 



24V27r 

> 1 — e 2 for all y. Clearly we obtain 



again from the first line. 



For the second assertion note that solving equation 



Co = ^ VT = Soml<y'''Vf) - <P{-la^^VT)) 
V 27r ^ ^ 



for given > yields the Black-Merton-Scholes implied volatility cr^^ . We obtain 
similarly as above 

.B 



< cr^^ - — = cr^"' - 

Do 



BS 



T 2 



BS^/f-^ 



^ 2 



T)) 



277^ 1 



< 



27r 
1 



(<l>(|)-'i>(-|)))L=.«. 



VT 



\/T 12^/2^ 



BSN3 



12 



□ 



Proposition [T] and [2] yield in particular the well-known asymptotic behaviour 
of an at the money call price in the Black-Merton-Scholes model for T — ^ cxd as 
described in 5]. 

Proposition [2] tells us that for the case when (o-s/T) <C 1 (which typically holds 
true in applications) , formula (|2.6|) gives a satisfactory approximation of the implied 
Black-Merton-Scholes volatility, and is very easy to calculate. We note that for the 
data reported by Bachelier (see [2] and [H]), the yearly relative volatility was of 
the order of 2.4% p. a. and T in the order of one month, i.e T = years so that 
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w 0.3. Consequently we get {aVTf w (0.008)^ w 5 x 10"'^. The estimate in 
Proposition [2] yields a right hand side of ^^'^^ 5 x 10"'' ~ 1.6 x lO^^S'o, i.e. the 
difference of the Bachelier and Black-Merton-Scholes price (when using the same 
volatility a — 2.4% p. a.) is of the order 10~^ of the price Sq of the underlying 
security. 

The above discussion pertains to the limiting behaviour, for aVT — > 0, of the 
Bachelier and Black-Merton-Scholes prices of at the money options, i.e. , where 
So = K. One might ask the same question for the case S'o 7^ K. Unfortunately, 
this question turns out not to be meaningful from a financial point of view. Indeed, 
if we fix 5*0 7^ if and let (7\/T tend to zero, then the Bachelier price as well as the 
Black-Scholes price tend to the pay-off function (S'o — of order higher than 
(aVT)", for every n > 1. Hence, in particular, their difference tends to zero faster 
than any power of iaVT). This does not seem to us an interesting result and 
corresponds - in financial terms - to the fact that, for fixed Sq K and letting 
a\/T tend to zero, the hinked pay-off function {Sq — K)^ essentially amounts to 
the same as the linear pay-off functions Sq — K, in the case Sq > K, and 0, in 
the case Sq — K < 0. In particular the functional dependence of the prices on 
a\/T is not analytical. This behaviour is essentially different from the behaviour 
at the money, S'o = K, since there the prices depend analytically on a^/T and the 
difference tends to zero of order 3. 



3. Further results of L. Bachelier 

We now proceed to a more detailed analysis of the option pricing formula (|2.2p for 
general strike prices K. Let C = Cq denote the option price from (|2.2p . We shall 
introduce some notation used by L. Bachelier for the following two reasons: firstly, 
it should make the task easier for the interested reader to look up the original texts 
by Bachelier; secondly, and more importantly, we shall see that his notation has 
its own merits and allows for intuitive and economically meaningful interpretations 
(as we have already seen for the normalization H = of the volatility, which 

equals the time-standardized price of an at the money option). 

L. Bachelier found it convenient to use a parallel shift of the coordinate system 
moving Sq to 0, so that the Gaussian distribution will be centered at 0. We write 

B Fp 

(3.1) a= P_ , m,:=K-So, P:=m + C. 

V 27r 

The parameter a equals, up to the normalizing factor -^=, the time standardized 

absolute volatility (j^^/T at maturity T. Readers familiar, e.g. with the Hull- White 
model of stochastic volatility, will realize that this is a very natural parametrization 
for an option with maturity T . 

In any case, the quantity a was a natural parametrization for L. Bachelier, as it 
is the price of the at the money option with the maturity T (see formula 12. 6p . so 
that it can be directly observed from market data. 

The quantity m is the difference between the strike price K and So and needs 
no further explanation. P has a natural interpretation (in Bachelier 's times it was 
called "ecart", i.e. the "spread" of an option): it is the price P of a european put 
with maturity T and strike price K, as was explicitly noted by Bachelier (using. 
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of course, different terminology). In today's terminology this amounts to the put- 
call parity. Bachelier interpreted P as the premium of an insurance against prices 
falling below K. 

This is nicely explained in [2]: a speculator "a la hausse", i.e. hoping for a 
rise of St may buy a forward contract with maturity T. Using the "fundamental 
principle" , which in this case boils down to elementary no arbitrage arguments, one 
concludes that the forward price must equal 5*0, so that the total gain or loss of 
this operation is given by the random variable St — Sq at time T. 

On the other hand , a more prudent speculator might want to limit the maximal 
loss by a quantity K > 0. She thus would buy a call option with price C, which 
would correspond to a strike price K (here we see very nicely the above mentioned 
fact that in Bachelier's times the strike price was considered as a function of the 
option price C ~ la "prime" in french - and not vice versa as today). Her total gain 
or loss would then be given by the random variable 

{St - K)+ - C. 

If at time T it indeed turns out that St > K , then the buyer of the forward 
contract is, of course, better off than the option buyer. The difference equals 

{St - So) -^[{St~K)~C]^K-So + C = P, 

which therefore may be interpreted as a "cost of insurance" . If St < K, we obtain 

{St - So) ~[0~C]^{St-~K) + K-Sq + C = {St - K) + P. 

By the Bachelier's "fundamental principle" we obtain 

P = E[{St - K)^]. 

Hence Bachelier was led by no-arbitrage considerations to the put-call parity. For 
further considerations we denote the put price in the Bachelier model at time t ~ 
by P(P or P{m) respectively. Clearly, the higher the potential loss C is, which the 
option buyer is ready to accept, the lower the costs of insurance P should be and 
vice versa, so that we expect a monotone dependence of these two quantities. 

In fact, Bachelier observed that the following pretty result holds true in his model 
(see [3], p.295): 

Propositions (Theorem of reciprocity) . For fixed > andT > the quantities 
C and P are reciprocal in Bachelier's model, i.e. there is a monotone, strictly 
decreasing and self-inverse (that is I = function I : IR>o R>o such that 

P = I{C) . 

Proof. Denote by ip the density of St — Sq, then 

/•oo 

C{m) = / {x — m)^p{x)dx, 

J ni 

/m 
{m — x)ip{x)dx. 
-oo 

Hence we obtain that C{—m) — P{m). We note in passing that this is only due to 
the symmetry of the density -0 with respect to reflection at 0. Since C'{m) < (see 
the proof of Proposition [T]) we obtain P — P{m{C)) := /(C), where C i— > m{G) 
inverts the function m ^ C{m). C maps R in a strictly decreasing way to M>o and 
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P maps M in a strictly increasing way to M>o. The resulting map / is therefore 
strictly decreasing, and - due to symmetry - we obtain 

I{P) = P{m{P)) = P{-m,{C)) = C, 

so / is self-inverse. □ 

Using the above notations (|3.ip . equation (|2.2p for the option price C(f (which 
we now write as C(m) to stress the dependence on the strike price) obtained from 
the fundamental principle becomes 

/•oo 

(3.2) C{m) = / {x-m)^JL{dx), 

where /i denotes the distribution of St — Sq, which has the Gaussian density fi{dx) = 
ip(x)dx, 

As mentioned above, Bachelier does not simply calculate the integral (|3.2p . 
He rather does something more interesting (see [31 p. 294]): "Si Ton developpe 
I'integrale en serie, on obtient ", i.e. "if one develops the integral into a series one 
obtains..." 

(3.3) C{m) = a \ —^ — + 

In the subsequent theorem we justify this step. It is worth noting that the method 
for developing this series expansion is not restricted to Bachelier's model, but holds 
true in general (provided that ^, the probability distribution of St — Sq, admits a 
density function tp, which is analytic in a neighborhood of 0). 

Theorem 1. Suppose that the law /i of the random variable St admits a density 

^{dx) = tl:{x)dx, 

such that tp is analytic in a ball of radius r > around 0, and that 

xip{x)dx < oo. 

J —oo 

Then the function 

/"OO 

C(m) — {x — m)^{dx) 

J m 

is analytic for |77i| < r and admits a power series expansion 

oo 

C(to) = y^Ckm'', 

where cq ~ x'4>{x)dx, ci — — ilj(x)dx and Ck — jftp^'^~^\0) for k >2. 

Proof. Due to our assumptions C is seen to be analytic as sum of two analytic 
functions, 

/"OO /"OO 

C{m) = / xil){x)dx — m / ip{x)dx. 
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Indeed, iftp is analytic around 0, then the functions x i~+ xip^x) and m i-^ x'ip{x)dx 
are analytic with the same radius of convergence r. The same holds true for the 
function m i— s- m 'ilj{x)dx. The derivatives can be calculated by the Leibniz rule, 



C' {m) = —mip{m) — / i/j{x)dx + 1711^(171) 
'ip{x)dx, 



C" irn) — V'(to), 
whence we obtain for the fc-th derivative, 

for fc > 2. □ 

Remark 1. IJ we assume thatm 1— *■ C{m) is locally analytic aroundm — (without 
any assumption on the density ip), then the density x 1-^ i^ix) is analytic around 
x ^ 0, too, by inversion of the above argument. 

Remark 2. In the case when -0 equals the Gaussian distribution, the calculation 
of the Taylor coefficients yields 

:(:^exp(-— ^)) - — ^ 



dy 2iTa Aira^ 4 tt^o'^ 

d^ ^ 1 ^^^^ y2 ^ 127ra2-?/2^_ 

dy"^ 27ra Ana^ 8 Tr^a^ 

— (— exp(-^^)) = 1 6TTa^y-y^ ^- 

dy^ 27ra Ana'^ 16 tt^o^ ' 

d'^ I y2 1 Utt'^o'^ -m/na'^ +y^ 

d7^2^''''P^~W^^ " 32 ^ 



1 » 



1 y 



1 !/ 



(7onsegMen% 0(0) = ^i^, V'(0) = 0, 0"(O) - -4;^^, 0"'(O) = and ^"{0) = 
f hence with C(0) = a and C"(0) = - 

^("^) = '^-y + 4^-96^ + T92(W+^^^^^ 
and t/ie series converges for allm, as the Gaussian distribution is an entire function. 
This is the expansion indicated by Bachelier in [3]. Since P{~m) — C{m), we also 
obtain the expansion for the put 

9 A fi 



(3.5) P{m) = a + \ ^ + ^ + 0(to«). 



TO 

y ' Ana 967r2a3 ' 192077%^ 

Remark 3. Looking once more at Bachelier 's series one notes that it is rather a 
Taylor expansion in ^ than in to. Note furthermore that — is a dimensionless 
quantity. The series then becomes 

m 

C{m)^aF{-) 
a 

We note as a curiosity that already in the second order term the number tt appears. 
Whence - if we believe in Bachelier 's formula for option pricing - we are able to 
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determine vr at least approximately (see 13.6)) below) - from financial market data 
( compare Georges Louis Leclerc Comte de Buffon's method to determine tt by using 
statistical experiments). 

Let us turn back to Bachelier's original calculations. He first truncated the 
Taylor series (|3.4p after the quadratic term, i.e. 

(3.6) + W 

This (approximate) formula can easily be inverted by solving a quadratic equation, 
thus yielding an explicit formula for to as a function of C. Bachelier observes 
that the approximation works well for small values of (the cases relevant for 
his practical applications) and gives some numerical estimates. We summarize the 
situation. 

Proposition 4 (Rule of Thumb 1). For given maturity T , strike K and > 0, 
let m = K — So and denote by a = C(0) the Bachelier price of the at the money 
option and by C{m) the Bachelier price of the call option with strike K — So + m. 
Define 

(3.7, eM.„-| + £l 

then we get an approximation of the Bachelier price C(m) of order 4, i.e. C{m) — 
dim) = 0{m^). 

Remark 4. Note that the value of C{m) only depends on the price a — C{0) of 
an at the money option (which is observable at the market) and the given quantity 
m = K — So- 

Remark 5. Given any stock price model under a risk neutral measure, the above 
approach of quadratic approximation can be applied if the density ip of St — So is 
locally analytic and admits first moments. The approximation then reads 

(3.9) C(to) = C(0)-B(0)TO+^m2 

up to a term of order 0(m^). Here C(0) denotes the price of the at the money 
european call option (pay-off {St — So)+ ), B{0) the price of the at the money 
binary option (pay-off i{ST>So}) '^'^'^ "0(0) the value of an at the money "Dirac" 
option (pay-off dso, with an appropriate interpretation as a limit). Notice that B{0) 
can also be interpreted, in Bachelier's model, as the hedging ration Delta. 

Example 1. Take for instance the Black- Merton-Scholes model, then the terms of 
the quadratic approximation iS. 9\) can be calculated easily, 

c(o) = Som^a^'Vr) - ci>(-ia^sVr)) 

B{0) = ^-^a^'VT), 
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Notice that the density ijj of St ~ >S'o in the Black-Merton-Scholes model is not an 
entire function, hence the Taylor expansion only converges with a finite radius of 
convergence. 

Although Bacheher had achieved with formula p.7p a practically satisfactory 
solution, which allowed to calculate (approximately) m as a function of C by only 
using pre-computer technology, he was not entirely satisfied. Following the refiexes 
of a true mathematician he tried to obtain better approximations (yielding still 
easily computable quantities) than simply truncating the Taylor scries after the 
quadratic term. He observed that, using the series expansion for the put option 
(see formula [23]) 

, m m? nri^ 
Pirn) = a H 1 ^ + ... 

and computing the product function C{m)P{m) or, somewhat more sophisticatedly, 
the triple product function C{m)P{m) ^ one obtains interesting cancel- 

lations in the corresponding Taylor series. 



C(m)F(m) =a^ - — + — + 0{m^), 
A 2-K 

^, sC(m)+P(m) n ni?a 'irn^a 

C{m)P{m) \ = « — + + ^("^ )• 

Observe that {C{m)P{m))^ is the geometric mean of the corresponding call and 
put price, while {C{m)P{m) )| jg tJ^g geometric mean of the call, the put 

and the arithmetic mean of the call and put price. 
The latter equation yields the approximate identity 

(C(to) + P{m))C{m)P{m) w 2a^ 

which Bachelier rephrases as a cooking book recipe (see |3!, p. 201): 

On additionne I'importance de la prime et son ecart. 
On multiplie I'importance de la prime par son ecart. 
On fait le produit des deux resultats. 
Ce produit doit etre le meme pour toutes Ics primes 
qui ont meme echeance, 

i.e. "One adds up the call and put price, one multiplies the call and the put price, 
one multiplies the two results. The product has to be the same for all premia, which 
correspond to options of the same maturity a." This recipe allows to approximately 
calculate for m 7^ m' in a quadruple 

{C{m),P{m),C{m'),P{m')) 

any one of these four quantities as an easy (from the point of view of pre-computer 
technology) function of the other three. Note that, in the case m = 0, we have 
C(0) = P(0) = a, which makes the resulting calculation even easier. 

We now interpret these equations in a more contemporary language (but, of 
course, only rephrasing Bachelier's insight in this way). 

Proposition 5 (Rule of Thumb 2). For given T > 0, > and m = K ~ So 

denote by C{m) and P(m) the prices of the corresponding call and put options in 
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the Bachelier model. Denote by 

a(m) := C{m)P{m) 

b{m) = C{m)P{m){ ) 

the products considered by Bachelier, then we have a(0) = and 6(0) — , and 

^.:-<^..,„.,, 
^-^^-<™'.- 

Remark 6. We rediscover an approximation of the reciprocity relation of Proposi- 
tion \3.1\ in the first of the two rules of the thumb. Notice also that this rule of thumb 
only holds up to order (■^)^. Finally note that tt — 3 ~ 0.1416 while tt — 2 « 1.1416, 
so that the coefficient of the quadratic term in the above expressions is smaller for 
by a factor of 8 as compared to ■ This is why Bachelier recommended this 
slightly more sophisticated product. 



4. Bachelier versus other models 

Bachelier's model yields very good approximation results with respect to the 
Black-Scholes-model for at the money options. This corresponds to notions of weak 
(absolute or relative) errors of approximation in short time asymptotics, i.e. esti- 
mates of the quantitiy 

|i?(/(5|))-i?(/(5|«))| 

for T — !■ 0. Compare |6j for the notion of weak and strond errors of approximation 
in the realm of numerical methods for SDE. This can be very delicate as seen in 
Proposition [2] and the remarks thereafter. However, in this section we concentrate 
on the easier notion of L^-strong errors, which means estimates on the L^-distance 
between models at given points in time T. 

We ask two questions: first, how to generalize Bachelier's model in order to ob- 
tain better approximations of the Black-Scholes-model and second, how to extend 
this approach beyond the Black-Scholes model. Both questions can be answered by 
the methods from chaos expansion. There are possible extensions of the Bachelier 
model in several directions, but extensions, which improve - in an optimal way ~ 
the L^-distance and the short-time asymptotics (with respect to a given model), 
are favorable from the point of view of applications. This observation in mind we 
aim for best (in the sense of L^-distance) approximations of a given general process 
{St)o<t<T by iterated Wiener-Ito integrals up to a certain order. We demand that 
the approximating processes are martingales to maintain no arbitrage properties. 
The methodology results into the one of chaos expansion or Stroock- Taylor The- 
orems (see [S]). Methods from chaos expansion for the (explicit) construction of 
price processes have already proved to be very useful, see for [7] . In the sequel we 
shall work on a probability space (fi, J^, P) carrying a one-dimensional Brownian 
motion (Wt)Q<-j<-j, with its natural filtration {^t)Q.^t<T- ^^"-"^ necessary details 
on Gaussian spaces, chaos expansion, n-th Wiener chaos, etc, see [7]. Certainly, the 
following considerations easily generalize to multi-dimensional Brownian motions. 

We shall call any Gaussian martingale in this setting a (general) Bachelier model. 
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Definition 1. Fix N > 1. Denote by (A'fj"'')o<t<T martingales witti continuous 
trajectories for < n < N , such tfiat M^"'^ G Tin for < t < T , where 7i„ denotes 
the n-th Wiener chaos in L^(ri). Then we call the process 

N 
ri=0 

an extension of degree N of the Bachelier model. Note that M^^^ — M*^°) is 
constant, and that S^^^ = M(o) + a4^^ is a (general) Bachelier model. 

Given an L^-martingale {St)o<t<T with Sq > and > 1. There exists a unique 
extension of degree N of the Bacheher model (in the sense that the martingales 
(M("'')o<t<T are uniquely defined for < n < N) minimizing the L^-norm E[{St — 
Sl^^y] for all < t < T. Furthermore S*}^^ ^ 5* in the L^.^orm as N ^ oo, 
uniformly for < t < T. 

Indeed, since the orthogonal projections p„ : L^{Q,!Ft, P) L'^{^,^t, P) onto 
the n-th Wiener chaos commute with conditional expectations E{.\!Ft) (see for 
instance |H|), we obtain that 

Af(") :=P„(5i), 

for < i < r, is a martingale with continuous trajectories, because 

E{pn{ST)\Tt) ^ pn{E{ST\Tt) ^ Pn{St) 
for Q <t <T. Consequently 

is a process minimizing the distance to St for any t G [0,T]. Clearly we have that 

oo 

in the L'^-topology. 

Example 2. For the Black- Merton-Scholes model with a^^ ^ a we obtain that 

\lt 

where Hn denotes the n-th Hermite polynomial, i.e. [n + l)i?„+i(x) — xHn(x) — 
Hn-i{x) and Hq{x) = 1, Hi{x) = x for n > 1 (for details see Hence 

We recover Bachelier 's model as extension of degree 1 minimizing the distance to 
the Black-Merton-Scholes model. Note that we have 
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for < t < T . Furthermore we can calculate a sharp constant Cn , namely 



c%= sup y: (-v^)^^— 

— — m>N-\-l 



^ ^ ' ml 

Tn>N+l 



0. 



since E[{H,n{^yf] ^ for < t < T and m > 

Due to the particular structure of the chaos decomposition we can prove the 
desired short-time asymptotics: 

Theorem 2. Given an -martingale {St)o<t<T, assume that St = Si^o^^(-^») 
with symmetric -functions fi : [0,T]' — > R and iterated Wiener-Ito integrals 



Wt\f^)■.^ / f^iti,...,t^)dWt,---dWt,. 

JO<ti<-<ti<t 

If there is an index io such that for i > io the functions fi are bounded and := 

Si>jo ^^\^ < ^^e*^ we obtain \\St - S't^-'||2 < Kt-^ for each N > io 

and 0<t<T. 

Proof We apply that E[iEiWUf^)\^tm = Wmf^\\hi[o.TY) < ItII/^II^ ^r 
i > io. Hence the result by applying 

\\St - 4'^'||2 < V E[{EiWi.{f,)\TT)f] < Kt"^. 



i>io 



□ 



Remark 7. Notice that the Stroock- Taylor Theorem (see [8 , p. 161) tells that for 
St G T>^'°" the series 

oo 

Y,W^{{ti,...,t,)^ E{Dt„...,uST)) = St 

i=0 

converges in T>^'°° . Hence the above condition is a statement about boundedness 
of higher Malliavin derivatives. The well-known case of the Black- Merton-Scholes 
model yields 

n qBS 1 (gii 

^ti,...,ti^T — -'-[o,T] 

for i > 0, so the condition of the previous theorem is satisfied. 
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5. Appendix 

We provide tables witli Baclielier and Black-Scholes implied volatilities for differ- 
ent times to maturity, where we have used Nasdaq-index-option data in discounted 
prices in order to compare the results. Circles represent the Black-Scholes implied 
volatility of data points above a certain level of trade volume. The dotted line 
represents the implied Bachelier (relative) volatility — Soarei- 



Black-Scholes and Bachelier (clotted) Volatility for an Index option maturing In 10 days 
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Black-Scholes and Bachelier (dotted) Volatility for an index option maturing in 73 days 
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